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Analysis of Energetic Models for 
Rate-Independent Materials 



Alexander Mielke* 

Abstract 

We consider rate-independent models which are defined via two function- 
als: the time-dependent energy-storage functional I : [0, T] x X — > [0, oo] and 
the dissipation distance 23 : X x X — ► [0, oo]. A function z : [0, T] — > X is 
called a solution of the energetic model, if for all < s < t < T we have 

stability: l(t,z(t)) < X{t,z) + T>(z(t),z) for all z £ X; 

energy inequality: T(t, z(t))+T)issT>(z, [s,t]) < 1(s, z(s))+ f* d T T(r, z(r)) dr . 

We provide an abstract framework for finding solutions of this problem. 
It involves time discretization where each incremental problem is a global 
minimization problem. We give applications in material modeling where z £ 
Z C X denotes the internal state of a body. The first application treats 
shape-memory alloys where z indicates the different crystallographic phases. 
The second application describes the delamination of bodies glued together 
where z is the proportion of still active glue along the contact zones. The third 
application treats finite-strain plasticity where z(t,x) lies in a Lie group. 
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1 Introduction 

Many evolution equations can be written in the abstract form 

Oe0¥(i(t)) + DZ(t,js(t)), (1.1) 

where z £ X is the state variable, I is the energy-storage functional, ^ : X — > [0, oo] 
is a convex dissipation functional, and d^S means the set-valued subdifferential (see 
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2 for this doubly nonlinear form). Rate-independency is realized by assuming that 
*]> is homogeneous of degree 1. 

We replace the above differential inclusion by a weaker energetic formulation, 
which is also more general since it allows for z-dependent dissipation functionals. 
For given I : [0J]xI-> [0, oo] and a given dissipation distance T> : X x X — ► 
[0, oo] satisfying the triangle inequality, we impose the energetic conditions of global 
stability (S) and the energy inequality (E) instead of (|l.lfl . A function z : [0, T] — > X 
is called a solution of the energetic model, if for all < s < t < T we have 

(S) T(t, z{t)) < X(t, z] + V(z(t), z) for all z G X; 
(E) T(t,z(t))+VissT>(z,[s,t])<l(s,z(s)) + J* d T I(r, z(t)) dr. 
Here, Dissp(z, [s, t]) is called the dissipation of z on the interval [s, t] and is defined 
as the supremum of Ylj=i ^(^fe-i)) z fe)) over all iV G N and all discretizations 
s = t < t\ < . . . < tjsr = t. 

Assuming D(zo,Zi) — ^(zi~zq), convexity of 1{t,-) and further technical 
assumptions, this energetic formulation is equivalent to (|1.1|) . see However, 
the latter form is more general as it applies to nonconvex problems and it doesn't 
need differentiability of 1 1— > z(t) nor of z <— > I(t, z). A related energetic approach to 
equations of the type ljl.l|l is presented in |20| . however, it remains unclear whether 
that method applies to the rate-independent case. 

In Section we discuss the abstract setting in more detail and in Section |21 
we provide existence results for solutions for given initial values z(0) = zq. The 
existence theory is based on time-incremental minimization problems of the form 

z k G argmin{ l(tk,z)+V(zk-i,z) \ z G X} 

and the BV bound for z : [0, T] — * X obtained via the dissipation functional satisfy- 
ing T>(zq, z\) > cd \\zq—zi\\. However, one needs additional compactness properties, 
if X is infinite dimensional. Here we propose a version where X satisfies coercivity 
with respect to an embedded Banach space Y, i.e., I(t,z) > — Ci+ci||z||y with 
ci, Ci, a > 0, where Y is compactly embedded in X. 

For the case of T> having the form T>{zq,z\) — ^(zi— zq) this theory was 
developed in ^B]. The case of general T> can be found in |10j. 

The flexibility of the energetic formulation allows for applications in continuum 
mechanics, where z : 57 — > Z plays the role of internal variables in the material 
occupying the body £1 C M d . Note that Z may be a manifold containing the 
internal variables like phase indicators, plastic or phase transformations, damage, 
polarization or magnetization. By Z we denote the set of all admissible internal 
states. The elastic deformation is (p : f2 — > M. d and T denotes the set of admissible 
deformations (p. 

Energy storage is characterized via the functional £ : [0, T] x JFx Z — > R, where 
t G [0,T] is the (quasi-static) process time, which drives the system via changing 
loads. In typical material models, £ has the form 

£(t,cp,z) = J Q W(x,T>cp(x), z(x))dx - (£ cxt (t) , ip) , 

where W is the stored-energy density and ^ext(^) denotes the external loadings. 
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Dissipation is characterized by an infinitesimal Finsler metric A : x TZ — > 
[0, oo], such that the curve z : [to, t\] —> Z dissipates the energy 

Diss (z, [to, ti]) — J n A(x, z(t, x), zit, x)) dxdt. 

The global dissipation distance T>(zo, z\) is then the infimum over all curves connect- 
ing zo with z\ . The relation to the abstract theory above is obtained by eliminating 
the elastic deformation via 

T{t, z) — inf{ E it, <f, z) | f G T } for z G Z and T(t, z) — +00 else. 

Obviously, the functional X is now fairly complicated and it is important to have 
rather general conditions in the abstract theory. 

In SectionQ]wc illustrate the usefulness of the abstract approach by discussing 
three quite different applications; however, the theory is used in other areas as well, 
e.g., in fracture mechanics @J[3| and in micro-magnetics [51 119|. 

Our first model describes phase transformations in shape-memory alloys as 
discussed in 1171 1181 Here z : £1 — > Z indicates either the microscopic 
distribution of the phases or a mesoscopic average of the microscopic distribution. 
In the first case we choose Z = Z p = {ei, . . . , e p } C MP, where e 3 denotes the j-th 
unit vector in M. p and in the second case we choose Z = conv Z p . In both cases the 
dissipation distance is given by a volume integral measuring the amount of volume 
which is transformed into another phase: T>[zo,Zi) — Jn A(zi(x)— Zq(x)) dx, where 
A : M p — > [0, 00 [ is convex and homogeneous of degree 1. This leads naturally to 
the basic space X = L 1 (f2,]R p ) and Z = {z e X \ z{x) S Z a.e. }. 

Including in £ an interfacial energy proportional to the area of the interfaces 
between regions of different phases provides a reduced energy X which is coercive 
in Y = BV(f2,K p ), see [5]. For an existence result in the case without interfacial 
energy we refer to |17| . 

The second application describes the delamination of a body SI which is glued 
together along n hypersurfaces Tj, j — l,...,n. The internal state z : T = 
L>iTj — > [0, 1] denotes the percentage of glue along T which remains in effect. 
The dissipation is given by a material constant ct> times the destroyed glue, i.e., 
T>(zo,zi) — ct> Jp zq(x)— z\(x) da(x) for z\ < zo and V(zo,z\) = +00 else. The 
basic underlying space is L 1 (r) and now compactness arises via the trace operator 
H 1 (ri) — > L 2 (r) which makes the reduced energy functional X weakly continuous. 

The final application is devoted to the modeling of elasto-plasticity with finite 
strains. There the internal variable z = (P,p) consists of the plastic transfor- 
mation P G SL(d) and hardening parameters p G K fc . Invariance under previ- 
ous plastic deformations leads to dissipation metrics which are left-invariant, i.e., 
A((P,p), (P,p)) — A((I,p), (P^ 1 P,p)). This geometric nonlinearity clearly shows 
that we need general dissipation distances T> avoiding any linear structure. In 
single-crystal plasticity A is piecewise linear in G s1(g0 which leads to Ba- 

nach manifolds and the dissipation metric is then a left-invariant Finsler metric. 
For applications in this context see 1131 IT2] . 
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2 Abstract setup of the problem 

We start with a Banach space X which is not assumed to be reflexive, since our 
applications in continuum mechanics (cf. Section 0J) naturally lead to spaces of the 
form L 1 (r2,K fc ). The first ingredient of the energetic formulation is the dissipation 
distance T> : X x X — > [0, oo] satisfying the triangle inequality: 

T^(zi,z 3 ) < V(zi,z 2 ) + T>(z 2 , z 3 ) for all z u z 2 , z 3 G X. 

We don't enforce symmetry, i.e., we allow for T>(zq, z±) ^ T>{z\,z§) as in Section 
14.21 We assume that there is a constant cp > such that T>(z , z\) > cx>\\zi— z a \\x 
for all zo,zi E X. The latter condition is in fact the one which determines the 
appropriate function space X for a specific application. Moreover, V is assumed to 
be s-weakly lower semicontinuous. (We continue to use the abbreviation s-weak for 
"sequentially weak".) We call T>(zq,zi) the dissipation distance from zq to z\. 
For a given curve z : [0, T] — > X we define the total dissipation on [s, t] via 

Dissp(z; [s, t}) = sup{ £f X>(«fa_i), z(t s )) | NeE, s=T <n< ■ • • <r N =t}. (2.1) 

The second ingredient is the energy-storage functional X : [0,T]xI^ [0, oo], 
which is assumed to be bounded from below and then normalized such that it takes 
only nonnegative values. Here t £ [0, T] plays the role of a (very slow) process time 
which changes the underlying system via changing loading conditions. For fixed time 
t, the map I{t 1 •) : X — » [0, oo] is assumed to be s-weakly lower semicontinuous, i.e., 
Zj — 1 z implies X(t , z) < lim infy^oo X(t, Zj ) . Moreover, we assume that for all z with 
2(t, z) < oo the function t 2(t, z) is Lipschitz continuous with \dtl(t, z)\ < Cj. 

Definition 2.1 A curve z : [0,T] — > X is called a solution of the rate-independent 
model (V,T), if global stability (S) and energy inequality (E) holds: 

(S) For all t G [0, T] and allzeX we have J(t, z(t)) < I(t, z) + V(z(t), z). 
(E) For all t , h with < t < ii < T we have 

J(ti,«(ti))+DfesD(*;[to,ti]) <I(t ,z(t )) + d t I(t,z(t))dt. 



The definition of solutions of (S)&(E) is such that it implies the two natural re- 
quirements for evolutionary problems, namely that restrictions and concatenations 
of solutions remain solutions. To be more precise, for any solution z : [0, T] — > E 
and any subinterval [s,t] C [0,T], the restriction z|r s ,t] solves (S)&(E) with initial 
datum z(s). Moreover, if z x : [0,t] E and z 2 : [t,T] -> E solve (S)fe(E) on the 
respective intervals and if Z\{t) = z 2 {t), then the concatenation z : [0,T] — ► E solves 
(S)&(E) as well. Under a few additional assumptions, it is shown in ^H] that (S) 
and (E) together imply that, in fact, the energy inequality is in an equality, i.e., for 
< t < ti < T we have 

J(ti,z(h))+Dvss v {z;[to,h])=X(to,z(to)) + d t T{t, z{t))dt. (2.2) 

Rate-independency manifests itself by the fact that the problem has no in- 
trinsic time scale. It is easy to show that 2 is a solution for (D,2) if and only 
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if the reparametrized curve z : t i— > z(a(t)), with a > 0, is a solution for (2?,Z), 
where T(t,z) = T{a{t), z). In particular, the stability (S) is a static concept and 
the energy estimate (E) is rate-independent, since the dissipation defined via (|2.1|l 
is scale invariant like the length of a curve. 

The major importance of the energetic formulation is that neither the given 
functionals T> and T(t, •) nor the solutions z : [0, T] — > X need to be differentiable. 
In particular, applications in continuum mechanics often have low smoothness. Of 
course, under additional smoothness assumptions on T> and T the weak energetic 
form (S)&(E) can be replaced by local formulations in the form of differential in- 
clusions like (f 1 . 1 1> ([21 ED]) or variational inequalities. See ^B] for a discussion of 
the implications between these different formulations. 

3 Time discretization and existence 

The major task is now to develop an existence theory for the initial value 
problem, i.e., to find a solution in the above sense which additionally satisfies z(0) = 
zq. In general, we should not expect uniqueness without imposing further conditions 
like smoothness and uniform convexity of T{t, •) and 2?, see [16| . 

The stability condition (S) can be rephrased by defining the stable sets 

S{t) :={zeX\ l(t, z) < l(t, z) + V(z, z) for all z e X }. 

Then, (S) simply means z(t) e S(t) for all t € [0, T]. The properties of the stable 
sets turn out to be crucial for deriving existence results. 

One of the standard methods to obtain solutions of nonlinear evolution equa- 
tions is that of approximation by time discretizations. To this end we choose discrete 
times = io<ti<...<iAr = T and seek z^ which approximates the solution z at 
tk, i.e., Zk ~ z(tk). Our energetic approach has the major advantage that the values 
Zk can be found incrementally via minimization problems. Since the methods of 
the calculus of variations are especially suited for applications in material modeling 
this will allow for a rich field of applications. 

To motivate the following incremental variational problem consider the non- 
linear parabolic problem h(dtu) — div(^4Du) + g, where we assume h'(v) > 0. The 
associated fully implicit incremental problem reads 

K t,-), A uk-Uk-i)) = div(ADw fc ) + g(t k ). 

v k v k— 1 

With H(v) = Jo h(w) dw we see that Uk must be a minimizer of the functional 
Jk(uk-i] ■):«!-> $ n {tk-tk~i)H( tk _ l tk _ 1 (u-Uk-i)) + f(ADu,Du) - g(t k )udx. 

In the simplest rate-independent case the function h is given by the signum function 
which implies H(v) = \v\. Hence, the length tk—tk-i of the fc-th time step disap- 
pears in the functional J^. In our more general setting the incremental problem 
takes the following form: 
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(IP) For zq £ X with 1(0, zq) < oo find z±, . . . , Zn € X such that 

Zk G argmmjX^fc, z) + V(zk-i, z) \ z £ X } for k = 1, . . . , N. (3.1) 



Here "argmin" denotes the set of all minimizers. Using the s-weak lower semi- 
continuity of D and X and the coercivity T(t, z)+T>(z k -i, z) > ct>\\z— Zk-i \\ we 
obtain the following result. 

Theorem 3.1 The incremental problem always has a solution. Each solution 
satisfies, for k = 1, . . . , N , the following properties: 

(i) Zk is stable for time t^, i-e., zj, 6 S(tk); 

(ii) f[ th _ lt t k ] d s X(s,z k )ds < l(t k ,z k ) - l(t k -i, Zk-i) +T>(z k -i,Zk) 

- f[t k .ut k ] dsX ^' Zk -^ ds; 
(in) T(t k , z k ) + £j =x V( Zj -i, Zj ) < 1(0, zq) + CjT; 
H \\z k \\ < \\z \\+(Z(0,zo)+CxT)/cv 

The assertions (i) and (ii) are the best replacements for the conditions (S) and (E) 
in the time-continuous case. 

For each discretization P = {0, t\, . . . , t^-i, T} of the interval [0, T] and each 
incremental solution (zk)k=i,...,N of (IP) we define two piecewise constant functions 
which attain the values Zk at tk and are constant in-between: Z p is continuous 
from the left and Z p is continuous from the right. Summing the estimates (ii) in 
Theorem 13. II over k = j, . . . ,m we find the following two-sided energy estimate. 

Corollary 3.2 Let P be any discretization of [0,T] and (zk)k=o,...,N a solution of 
(IP), then for < j < m < N we have the two-sided energy inequality 

l{t 3 ,Z p {tj)) + d s l{s, Z p {s))ds <l{t m , Z P (t m )) + T)iss v (Z P , [tj,t m ]) 

<T(t j ,Z p (t j ))+f*™ d s l(s,Zp(s))ds. 

The existence of solutions can now be established by taking a sequence (P(Z))/ e pj 
of discretizations whose fineness 6^ = max{ tj—t^^ | j = 1, . . . , } tends to 
0. Moreover we assume that the sequence is hierarchical with P(l) C P(l+1). The 
associated solutions of (IP)C' define z^ 1 ' := Z p ( l \ The construction of a solution of 
(S)&(E) consists now of two parts. 

First we use the dissipation bound (iii) of Theorem 13.11 to obtain an a priori 
bound in BV([0,T], X): 

c v J [Q T] \\dzW\\ x <Diss v (z^,[0,T])<l(0,z ) + C x T. 

Then, Part (iv) in Theorem l3.1l and the following additional compactness condition 
H3.2fl allows us to apply a selection principle in the spirit of Helly. 

For all R > and all t e [0, T] the sets 

TVp := { z S X | T>(zq, z) < R, X(t, z) < R} are s-weakly compact. 
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Thus, we can extract a subsequence (Z„)„ g n such that for all t G [0, T] the se- 
quence z^(t), n G N, converges weakly to a limit z^°°'(t) with Dissx>(z(°°\ [0, T]) < 
liminf„^ oo Diss I ,(z('"),[0,r]). 

Second we need to show that z^°°> is a solution of (S)&(E). Using Corollarv l3.2l 
it is easy to give conditions which guarantee that satisfies (E) for t = and 
t\ = T, and by (12. 2[) this is sufficient. To obtain stability of z^ 00 ^ there are essentially 
two different ways. If additional compactness properties allow us to conclude that 
the convergence of z^ ln ^(t) to z^ 00 ^ also happens in the strong topology, then we are 
in the good case. Then it suffices to know that the set 

S[o,t] = { (t, z) G [0, T] x X | z G S(t) } = U te[0 , T ] (t, S(t)) 

is closed in the strong topology. If strong convergence cannot be deduced, one needs 
to show that <5>[o,t] is s-weakly closed. This property is quite hard to obtain, since 
even under nice convexity assumptions on T(t, •) the sets S(t) are generally not 
convex. 

The following theorem provides two alternative sets of assumptions which en- 
ables us to turn the above construction into a rigorous existence proof. 

Theorem 3.3 Let T> and X be given as above and satisfy QtJ.ty) . If one of the 

conditions (a) or (b) is satisfied, then for each zq G X with 1(0, zq) < oo there is 
at least one solution z G BV([0,T],X) of (S) k (E) with z(0) = z . 

(a) The set <S[o,r] * s s-weakly closed and z i— ► dt1(t,z) is s-weakly continuous. 

(b) The sets TZ l R in are compact, the set <S[o,t] * s dosed, and z i— ► dp[(t,z) is 
continuous (all in the norm topology of X). 

Simple nontrivial applications of this theorem with either condition (a) or (b) 
are as follows: Let X = L 1 (f2) with C t d bounded and choose the dissipation 
distance T>(zq, z%) = ct>\\z\—zq\\x = cp J \z\(x)—Zo(x)\dx. As a first case consider 

T\(t,z) — J n a(x)\z(x)\P— g(t, x)z(x)dx + 7, 

where a(x) > a > 0, j3 > 1, and g G C 1 ([0, T], L°°(f2)). The sets are closed 
convex sets which lie in the intersection of an L 1 -ball and an Lr-ball. Hence, we 
obtain the s-weak compactness condition (|3.2() . Yet, IZr is not strongly compact in 
L 1 (fi). The stable sets for T\ are given by 

5 x (t) = {z G Ll(fi) I \z(x)\^z(x) G [^g^, ^gg^] for a.a. id!!}, 

which shows that they are s-weakly closed since they are convex and closed. Hence, 
condition (a) is satisfied. 

As a second case consider the nonconvex energy functional 

l 2 (t,z)= j n T;\T)z(x)\ 2 +f(t, x, z{x))dx for z G H 1 (il) and + 00 else, 

where / : [0, T] x £1 x M. — > R and dtf are continuous and bounded. Now, TZ^ is 
already compact in L 1 (J7) since it is closed and contained in an H 1 -ball. With these 
properties, it can be shown that condition (b) of Theorem 13 . 31 holds . 
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4 Applications in continuum mechanics 

The flexibility of the energetic formulation allows for applications in continuum 
mechanics. We consider an elastic body which is given through a bounded domain 
Slcl 1 ' with sufficiently smooth boundary. The elastic deformation is given by the 
mapping tp : ft — > M. d , and the set of all admissible deformations is denoted by T, 
which implements the displacement boundary conditions. 

The variable z e Z includes all the internal variables like phase indicators, 
plastic or phase transformations, damage, polarization or magnetization. A function 
z : ft — > Z gives the internal state of the material, and Z denotes the set of 
all admissible internal states. Note that Z may be a manifold with (nonsmooth) 
boundary. In plasticity we have Z = SL(rf) x R fe , in phase transformations we let 
Z = {z <E [0, l] fe | Y7i z ^ = 1 }i an d m micro-magnetism z is the magnetization 
satisfying \z(t, x)\ = m > 0. Moreover, below we will also consider an application 
where z is not defined on all of ft but at certain parts of the boundary. 

Energy storage is characterized via the functional £ : [0,T]xf xZ which 
is the sum of the total elastic energy and the potential energies due to exterior 
loadings (Gibbs' energy): 

£{t, <f>, z) = f n W{x, Dcf>{x), z(x))dx - (4xt(i), <f>). 

Here t € [0, T] is the (quasi-static) process time which drives the system and the 
external loads are (£ ext (t),(f>) = J n f ext (t, x) ■ (f>(x)dx + / rtract g cx t(t, x) ■ (j>(x)da(x). 

Dissipation is characterized via the metric A : ft x TZ — > [0, oo] such that the 
curve z : [0, T] — > Z dissipates the energy 

Diss(z, [t ,*i]) = Jj* 1 J n A(x,z(t,x),z(t,x))dxdt on[t 1 ,t 2 \. 

For each material point x € ft, the infinitesimal metric A(a;, ■,■) : TZ — > [0, oo] 
defines a global distance function D(x, •, •) : Z x Z — ► [0, oo] and on Z we obtain 
the global dissipation distance 

V(z ,zi)= J n D(x,z (x),Zi(x))dx 

= inf{ Diss (z, [t , h])\z€ C Li P([0, l],Z), z(0) = z , z(l) = z x }. 

The rate- independent problem for this material model is defined as in the above 
abstract part, but now the elastic deformation appears as an additional variable, 
which, however, does not generate any dissipation. 

Definition 4.1 A pair (<p,z) : [0, T] — > T x Z is called a solution of the rate- 
independent problem associated with V and £ if the global stability (S) and the 
energy inequality (E) hold: 

(S) For all t e [0, T] and all ($,z)eJ 7 xZ we have 

£(t,(f>(t),z(t)) < £{t,4>,z)+V{z{t),z). 
(E) For all t , h with < t < h < T we have 

£{t u 0(ti), z(ii))+Dis Sl7 (z; [t , h}) < £(t , <f>(t ),z(t ))+ t d t £(t, cf>(t), z(t))dt. 
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The connection with the above abstract theory is obtained by minimization 
with respect to the deformations 4> € J-, since the stability condition implies that 
4>{t) must be a minimizer of £(t, z(t)). We define the associated X via 

I(t, z) = inf { £ (£, ip, z) | tp € J 7 } for zeZ and + oo else. 

While this elimination is suitable for an abstract treatment, the practical approx- 
imation of solutions via the incremental approach is better done by keeping the 
deformation and eliminating the internal variable in each incremental step. In fact, 
in (IP) we now have to find 

O fe ,z fc ) G axgmin{£(t k ,'$,z) +T>(z k -u z) \ ($,z) £ T x Z }. (4.1) 

In this minimization problem the internal variable occurs only locally under the 
integral over f2 and hence can be eliminated pointwise. Defining the local reduced 
constitutive functions 

* rcd (z old ; x, F) := min{ W(x, F, z) + D(x, z old , z)\z£Z}, 
Z n cw(z old ; x, F) e argmin{ W(x, F, z) + D(x, z u, z) \ z £ Z}, 

and the reduced functional £ red (z i<i; t, (p) = J n & (z \<i; D<j>) dx — (£ C xt{t), 4 1 ) t ne 
solution of (|4.1|l is equivalent to finding <j> £ argminj £ Tcd (z k _i; t k , <j>) \ t/> £ T} and 
then letting z k = Z ncw (zfe-i; D0 fc ). For more details we refer to [T^. 

4.1 Phase transformations in shape-memory alloys 

We assume that, in each microscopic point y, an elastic material is free to 
choose one of p crystallographic phases and that the elastic energy density W is 
then given by Wj(D<fi). If the model is made on the mesoscopic level, then the 
internal variables are phase portions z^' £ [0, 1] for the j-th phase. We set Z = 
{z £ [0, If C R p | YX z(j) = 1} and X = L 1 (f2,R p ). The material properties 
are described by a mixture function W : R dxd x Z — > [0, oo], see j^EIEl- The 
dissipation can be shown to have the form D(zq,zi) = ip(zi-zo) with ip( v ) — 
max{ a m ■ v \ m — 1,...,M} > C^|w|, where a rn £ M. p are thermodynamically 
conjugated threshold values. 

So far we are unable to prove existence results for this model in its full gener- 
ality. However, the case with only two phases (p = 2) has been treated in [TJ under 
the additional assumption that the elastic behavior is linear and both phases have 
the same elastic tensor. In that case, one sets z = (8, 1—9) with 9 £ [0, 1]. It can be 
shown that I is a quadratic functional in 9 £ L 1 (f2, [0, 1]) C L 2 (f2). It then follows 
that the compactness condition (|3.2|l holds and condition (a) in Theorem 13.31 can 
be verified using the H-measure to handle the weak convergence of the nonconvex 
terms. 

A microscopic model is treated in There no phase mixtures are allowed, 
i.e., we assume z £ Z p := {e\,e^ ■ ■ ■ ,e p } C M p , where ej is the j-th unit vector. 
Thus, the functions z £ Z are like characteristic functions which indicate exactly 
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one phase at each material point. The dissipation is assumed as above, but now 
the elastic energy contains an additional term measuring the surface area of the 
interfaces between the different regions: 

£(t,<f>,z) = f n W(Dcl>,z)dx + af n \Dz\-(£ ext (t),<l)), 

where a is a positive constant and f Q \Dz\ is \/2 times the area of all interfaces. 
Here Z = { z : Q — > Z p \ J n \Dz\ < oo } and we set £(t 7 (j), z) = +00 for z £ Z. 

Hence, after minimization with respect to <fi we still haveZ(i, z) > 7+er J n \Dz\. 
This term provides for V} R (cf. an a priori bound in BV(f2,M p ) and hence 

we conclude compactness in X — L 1 (il,R p ). Under the usual additional con- 
ditions for the elastic stored-energy densities Wj we obtain for each zo 6 Z a 
solution (0,z) with (f> G L^(]0,T[,W 1 ' 2 (O,M <i )) and z € BV([0, T], L 1 ^, W)) n 
L^,(]0,T[,BV(r2,R p )) with z(t) € Z for all i € [0,T], see 0. 

4.2 A delamination problem 

Here we give a simple model for rate-independent delamination and refer to 
for a better model and the detailed analysis. 

Consider a body f2 C R d which is given by an open, bounded, and path- 
connected domain. Assume that the interior of the closure of differs from fi by 
a finite set of sufficiently smooth hypersurfaces Tj, j — 1, ... ,71. This means that 
with T := Uj=i we have int(cl(0)) = fiUT. We assume that the two sides of the 
body are glued together along these surfaces and that the glue is softer than the 
material itself. Upon loading, some parts of the glue may break and thus lose its 
effectiveness. The remaining fraction of the glue which is still effective is denoted 
by the internal state function z : T — > [0,1]. 

We let Z = { z : T —> [0,l]\z measurable} C X = L 1 ^). The dissipation 
distance T>(zq, z%) is proportional to the amount of glue that is broken from state 
zo to state z\\ 

T>(zo,Zi) — cd J r zo(y)—zi(y)da(y) for zq > Z\ and +00 else. 

Here we explicitly forbid the healing of the glue by setting V equal 00, if zq ^ z\. 

The energy is given by the elastic energy in the body, the elastic energy in the 
glue, and the potential of the external loadings: 

£(t,^z) = f n W(D<P)dx + f r z(y)Q(y, [0] r (l/))do(y) - (4* 

where for y € T the vector [0]r(y) denotes the jump of the deformation <f> across 
the interface T and Q(y, •) is the potential for the elastic properties of the glue. 

For simplicity we assume further that W provides linearized elasticity and Q 
is quadratic is well, then there is a unique minimizer <f> — <f>(t,z) € H 1 (r2,R d ) of 
£ (t, •, z). It can be shown that the mapping ■) : Z C L 1 (r) — > H x (0,R d ) is com- 
pact, which implies that the functional T(t, ■) : Z — > [0, oo[ is s- weakly continuous 
with respect to the L 1 -topology on Z. For the latter argument it is essential that 
z appears only linearly in the definition of £(t, <p, z). Theorem (I3.3f) with condition 
(a) provides the existence of solutions. 
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4.3 Elasto-plasticity 

The above theory can be applied to linearized elasto-plasticity, see Here 
we want to report on recent results concerning elasto-plasticity with finite strain. 
However, for this application the abstract existence theory is not yet available. 

Elasto-plasticity with finite strains is based on the multiplicative decomposi- 
tion of the deformation gradient F = in the form D0 = Fciast-P 1 where the 
plastic transformation P lies in the Lie group SL(d) = { P € R dxd | dot P = 1 }. The 
internal variable has the form z — (P,p) £ Z where pel' denotes the hardening 
parameters. We refer to P3 1181 IT2*] for mechanical motivations and mathematical 
details. For simplicity, we mention here only the case without hardening where 
z = P G SL(tf) =: Z and refer to f° r more general cases. 

The important point in finite-strain elasto-plasticity is that the dissipation dis- 
tance must be invariant under previous plastic deformations, i.e., D(QPq,QPi) = 
D(P , Pi) for all Q £ SL(ef). Equivalently, the infinitesimal metric A : TZ — > [0, oo] 
is left-invariant, i.e., A(P,P) — A(I, P~ 1 P). This implies that the dissipation dis- 
tance is characterized by a norm A(J, •) on sl(d) = TjSL(d) and that D(P ,Pi) 
behaves logarithmically in Pj which introduces strong geometric nonconvexi- 
ties. So far, even the solution of the incremental problem (IP) is not understood 
completely. Even in simple cases one has to expect non-attainment in (IP), which 
leads to the formation of microstructure. The easiest way to see the problems is to 
study the reduced energy density \l/ rcd in (|4.2(l . If this density is not quasi-convex, 
then there are loadings such that (IP) has no solution and relaxation techniques 
have to be employed, cf. |14|. 
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